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We summarize our recent theoretical studies on the spin dynamics in the nor- 
mal state of the metallic cuprates. The contrasting wave vector dependence of the 
dynamical spin structure factor S{q, uj) in LaSrCuO and YBaCuO systems are at- 
tributed to the differences in the fermiology, in conjunction with strong Coulomb 
correlations. These effects are found to account also for the anomalous temperature 
and frequency dependence of 5(q, uj). We conclude that the low energy spin dynam- 
ics of the metallic cuprates are described in terms of correlated quasiparticles with a 
Luttinger Fermi surface and a non-zero antiferromagnetic exchange interaction. 
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I. INTRODUCTION 



In the past few years, the study of various physical properties have provided many 
insights into the nature of both the normal and the superconducting states of high T c 
copper oxides. [^,0 We discuss here the spin dynamics in the metallic cuprates, with 
a focus on the normal state. The systematic study of the spin dynamics is motivated in 
part by the widespread belief that they may be relevant to the superconductivity. It is also 
of importance in elucidating the nature of doped Mott insulators in general, a fundamental 
and long standing problem in condensed matter physics. 

The spin dynamics in the cuprates were studied first in the undoped compounds, such 
as La2Cu04 and YBa2Cu306- The electronic structure in these systems is characterized by 
the half-filling of the outermost (copper) band. According to the one electron picture, these 
half-filled cuprates should be metals. However, they all show three-dimensional magnetic 
ordering at low temperatures, and continue to be insulators above the Neel temperature. 
In this paramagnetic insulating phase, spin dynamics have been studied at low energies 
through inelastic neutron scattering, and at higher energies from both (two-magnon) Raman 
scattering || and (zone-boundary one-magnon) neutron-scattering [[?[]. These results have 
led to a fairly convincing picture that, || the paramagnetic phase is described by the two- 
dimensional spin | Heisenberg model defined on the copper square lattice within a CuC>2 
layer. Combined with the existence of a large charge gap, || this indicates that the half-filled 
cuprates are Mott (Charge-transfer) insulators: the electrons are localized due to the strong 
on-site Coulomb repulsions, and exhibit at low energies spin degrees of freedom only. 

As the system is doped away from half filling, the magnetic ordering quickly vanishes, 
and an insulator to metal (superconductor) transition sets in. Our understanding of the half- 
filled case does not, however, uniquely specify the nature of the doped metallic phase. Many 
questions can be asked about the role of doping. For example, 1) doped holes disturb the spin 
background. Do they lead to a frustrated spin system, or manage to change the nature of 
the spins altogether? The observation of a Luttinger Fermi surface in the metallic cuprates 



2 



implies that, the localized copper spins of the half-filled limit are converted to itinerant 
electrons at low energies in the doped regime. How is this manifested in the spin dynamics? 
2) the Mott insulating nature of the half-filled case suggests that, the cuprates have strong 
Coulomb correlations. Is this relevant to the understanding of the low energy spin dynamics 
in the doped metallic phase? 3) in the half-filled case, the magnetic interactions are mainly 
between the nearest neighbor copper spins. What is the nature and the magnitude of 
the effective exchange coupling between low energy spin excitations in the doped metallic 
phase? Determining this coupling is of course important for the purposes of addressing 
the relevance of spin fluctuations to superconductivity; and 4) in the Mott insulating state, 
charge excitations are pushed to high energies while spin excitations remain at low energies. 
How do the doped holes modify the coupling between charge and spin degrees of freedom? 

In the metallic cuprates, a host of anomalies in the low energy spin dynamics have been 
revealed in nuclear magnetic/quadrupolar resonance (NMR/NQR) and inelastic neutron 
scattering experiments. One striking feature is the strong contrast between the spin fluctu- 
ation spectra, 5(q, u), in LaSrCuO and YBaCuO families. Other anomalies are observed 
in the wave vector, temperature and frequency dependence of 5"(q, u). In the following, we 
will demonstrate that the low energy particle-hole spin excitations near the Luttinger Fermi 
surface can naturally explain the differences in >S(q, uo) between LaSrCuO and YBaCuO 
systems. At the same time, strong Coulomb correlations play an important role for un- 
derstanding these magnetic data. Our analysis also establishes the existence of a non-zero 
antiferromagnetic exchange coupling between the quasiparticles. Finally, we will briefly 
comment on the relation between low energy spin and charge dynamics. 

The rest of the paper is organized as follows. In section |TJ we give a brief overview of the 
spin dynamical data in the metallic cuprates. Section [TTI| presents a formalism for magnetic 
interactions and dynamical spin susceptibility within the extended Hubbard model in the 
limit of strong Coulomb correlations. This formalism is used to analyze the spin fluctuation 
spectra in Section [TV]. Section |V] contains some brief discussion on the spin dynamics in 
the superconducting state and the relation between spin and charge dynamics. Several 
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concluding remarks are made in Section [V]]. 



This brief review summarizes mainly our own work, |]H , |1"2" ,1"3|,1"4"| and is not intended 
to be comprehensive. Readers are directed to other reviews for more discussions on the 
subject. [[| A more complete analysis of the normal state properties in the metallic cuprates 



along similar lines can be found in Ref. [ID], which also contains more extensive references. 
The importance of fermiology in understanding the spin fluctuation spectra in LaSrCuO 
has been emphasized independently by Littlewood et al. [15|] in which correlation effects are 
treated through inelastic quasiparticle lifetimes within the Marginal Fermi liquid scheme. 
The role of fermiology in understanding the contrasting spin dynamics in LaSrCuO and 
YBaCuO has also been emphasized subsequently by Tanamoto et al. |L6j within the gauge 
theory approach to the t— J model. These studies have all led to the conclusion that the bulk 
of the low energy spin fluctuation spectra are described in terms of itinerant renormalized 
quasiparticles, with no additional exotic dynamics in the spin channel. 



II. MAGNETIC DATA IN THE METALLIC CUPRATES 

Both NMR/NQR spin-lattice relaxation rate and the inelastic neutron scattering cross 
section measure the dynamical spin structure factor S(q,co). From the fluctuation- 
dissipation theorem, 

3(q,u) = (l+n&(w))Imx(q,w) (1) 

where x(q, u;) is the transverse dynamical spin susceptibility, and n&(a>) = l/(e u ^ T — 1) is a 
Bose factor. 

The spin- lattice relaxation rate, jr, is the rate at which the nuclear magnetization relaxes 
towards equilibrium starting from a non-equilibrium distribution. Such a spin flipping re- 
laxation is induced by the electronic spin excitations, through a hyperfine coupling between 
the nucleus and the electrons. Specifically, 

d)r ~ £^(q)s(q>^-o) (2 ) 
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Since the nuclear resonance frequency u is essentially zero in the scale of electronic en- 
ergies, jr measures spin dynamics at the vanishing frequencies. Due to the local nature 
of the nucleus dynamics, measures the wave vector q-integrated S(q,uj ), weighted by 
the q-dependent hyperfme coupling constant A(q). Qualitative information about the q- 
dependence in S(q,uj ) can be derived if there exist several kinds of nuclei at different sites 
(labeled by r in Eq. (^)) which relax via the same dynamical structure factor, providing 
the associated hyperfme coupling A r (q) has dominant contributions from different regions 
in the q-space for different sites r. 

Inelastic neutron scattering experiments measure the dynamical spin structure factor 
more directly. Here the cross section 

d 2 a 
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S(q,u) (3) 



The fact that both the wave vector and the frequency are resolved is of course an important 
advantage of these measurements as compared with NMR/NQR experiments. 

That the spin dynamics in the metallic cuprates are anomalous compared to conventional 
metals has been discussed extensively from the temperature dependence of To appre- 
ciate this anomalous behavior, we note that within a conventional metal can be loosely 
estimated as follows. The number of electrons available to flip a nuclear spin is proportional 
to TN(Ep), where N(Ep) is the density of states at the Fermi level. Each electron has 
a contribution to the relaxation rate proportional to N(Ep), which measures the number 
of available final states. Therefore, jr has a linear temperature dependence. Furthermore, 
since the static uniform spin susceptibility, Xi is proportional to the density of states N(E F ), 
there exists the Korringa relation, 

' 1 (4) 



T,T X 2 

(where dimensionless units are used.) 

Shown in Fig. [I] are the relaxation rates at both the planar copper and oxygen sites in 
(nearly optimal doped) YBa2CuOy_5 and that for the planar copper site in L^Sro.isCuCU. 



fi7| , |T8| j4|| Two peculiar features are associated with (jr)cu'- a) it is strongly enhanced com- 
pared to that in conventional metals; The nominal Korringa ratio is at least one order of 
magnitude larger; |TJ| and b) its temperature dependence deviates considerably from linear 
behavior. Finally, comparing (^j-)o and (jr)cu gives yet another important feature: c) (jr)o 
is almost linear in temperature, and has a nominal Korringa ratio of order 1. 

Within the "standard" one-spin-component interpretation, these anomalies are under- 
stood as follows. |]i^[l9ip0ipll , [l2| In a CuC>2 plane of the metallic cuprates, the spin degrees 



of freedom are mainly associated with the planar copper sites. Strong antiferromagnetic 
spin correlations are responsible for the excess relaxation of the copper nuclear magnetiza- 
tion, leading to an enhanced {^r)cu- The amount of enhancement depends on temperature. 
Therefore, {^r)cu deviates from the Korringa linear in T behavior. On the other hand, since 
the oxygen sites are located at the mid-point between two copper sites, the contributions to 
(fr) from neighboring copper spins cancel with each other. The anomalies are therefore 
absent for (^j-)o- 

In momentum space, one can infer that this interpretation is based on the assumption 
that the dynamical spin susceptibility is strongly peaked at the antiferromagnetic wave 
vector, Qaf = , -)■ The existence or not of such sharp q-structure in the dynamical spin 
susceptibility can, of course, be more directly established from inelastic neutron scattering 
experiments at low energies. In the past, neutron scattering studies have elucidated the 
nature of spin fluctuations in other strongly correlated fermion systems such as liquid 3 He 



and heavy fermion materials. p2 



The momentum dependence in the inelastic neutron scattering cross sections in metal- 



lic YBaCuO and in LaSrCuO are shown in Fig. |2|. p3| , p4|p^ , p6f The main features can 
be summarized as follows: a) in YBaCuO, antiferromagnetic spin fluctuations are essen- 
tially commensurate. However, the peaks are rather broad, and the peak widths are nearly 
temperature independent; b) in LaSrCuO, spin fluctuations are sharply peaked at incom- 
mensurate wave vectors Q*. The incommensurability, defined as the deviation of Q* from 
Qaf, quickly increases as the doping concentration is increased as shown in Fig. |](c). This 



incommensurability is much larger (by more than a factor of two) than that expected within 
an itinerant description in terms of the one band Hubbard model with nearest neighbor 
hopping only. 

In addition to the q-structure, the spin dynamics have also been extensively probed 
as a function of temperature and frequency. Within the lightly-doped non-superconducting 



LaSrCuO system |27| u ) is found to scale with ~, which is consistent with the Marginal 
Fermi liquid ansatz. f28j Systematic studies of the temperature and frequency dependence 
of x"(q, cu) at higher doping concentrations, while still under way for the normal state of 
superconducting LaSrCuO |p9 |, have however revealed low energy scales in YBaCuO as is 



illustrated in Fig. [3] and discussed extensively in Refs. p5| . [26 

It is indeed intriguing that, the two families of cuprates show contrasting q-dependence 
in the dynamical spin susceptibility. In the following, we show how this difference helps 
understand the nature of the low lying spin excitations. We will also demonstrate that the 
dynamical spin susceptibility revealed in inelastic neutron scattering in both YBaCuO and 
LaSrCuO can not be easily reconciled with the form- factor cancellation argument for (^r)o- 

III. STRONGLY CORRELATED QUASIPARTICLE DESCRIPTION 
OF THE MAGNETIC DYNAMICS 

To understand the spin dynamical spectra, we need a microscopic formalism for the 
dynamical spin susceptibility x(q, tu). In the following, we present such a scheme within the 



extended Hubbard model in the presence of strong Coulomb correlations. [ I I 



A. The Extended Hubbard Model in the Large U Limit 

We consider the extended Hubbard model defined in a Cu02 layer, 
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+ E Vpd{dlPi* + h.c.) +Y. U 4^4i d H (5) 

i I a i 

Here p and d denote the oxygen p x , p y and the copper d x 2_ y 2 orbitals respectively. We 
include a hybridization V p d between nearest neighbor copper and oxygen orbitals, a finite 
oxygen dispersion derived from oxygen-oxygen hopping matrix element t pp (nearest neighbor) 
and t' (next nearest neighbor). The oxygen and "bare" copper levels are called e p and e° d 
respectively. The on-site Coulomb repulsion for the copper is U. The parameters involved 
in the Hamiltonian (|^) appropriate to the copper oxides have been extensively studied. 



Because of the large number of parameters involved, we presume that the derived parameters 
only specify an appropriate range. Specific parameters within this range will be chosen from 
phenomenological constraints to be discussed shortly. Before doing this, we first give a 
qualitative discussion for the case that the Coulomb repulsion U and the level separation 
e p — e° d are large. 

For the half-filled case, the system is a Mott (Charge-transfer) insulator. Because of 
the exclusion of the double occupancy, the magnetic interaction has no static on-site (q- 
independent) component. The dominant contribution comes from the nearest neighbor 
term. The low energy effective Hamiltonian is the spin | Heisenberg model. Formally, 

■J(q) = Jo{cos(q x ) + cos(q y )) (6) 

where the nearest neighbor superexchange interaction J Q can be derived from the full model 
(§) through integrating out high energy dynamics. The relevant high energy states is mainly 
oxygen in character. This nearest-neighbor form is established to a good accuracy through 
high energy spin wave measurement. [[/]] 

In principle, there are many possible phases to which a Mott insulator may evolve when 
doped. The study of these different fixed points and the associated quantum phase tran- 
sitions is a fundamental and unsolved problem, which is beyond the scope of the present 
paper. We consider here the scenario that the electron localization and the antiferromagnetic 
ordering in the half-filled Mott insulator evolve into two corresponding incipient instabilities 
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when doped. |10 The incipient electron localization is characterized by an enhanced quasi- 
particle mass, while the incipient magnetic instability is manifested as an effective exchange 
coupling between the quasiparticles. Phenomenologically, such a picture is consistent with 
the existence of a Luttinger Fermi surface, the small and doping dependent plasma frequen- 
cies, and the presence of doping dependent magnetic fluctuations which we will analyze in 
detail. Further support for this picture comes from a comparison with anomalous physical 
properties in heavy fermions. |i0|,|3ll This picture is also consistent with indications from the 
small cluster diagonalization studies that, doping leads to the formation of quasiparticles 



with enhanced mass and a Luttinger Fermi surface. |32] We emphasize that, our results are 
more general than the particular approach we use to describe the model Hamiltonian Eq. 
(H). For our purposes here, this scheme can be viewed as providing a handle to differentiate 
itinerant, and strongly Coulomb renormalized, quasiparticle description of the spin dynamics 
in the metallic cuprates, versus other approaches. 

Formally, both the incipient localization p3| , |34] , |35|| and the magnetic interaction effects 
|TT| can be systematically studied within the large N approach to the model Hamiltonian ([S]) 
(where N is the spin degeneracy). This approach has been extensively used in the description 
of the Kondo lattice. For our purposes, the large N approach provides a convenient tool 
to systematically study the doping dependence of magnetic properties. In the end, it would 
be useful to compare our conclusions with numerical analysis in a small cluster on generalized 



Hubbard models with extended dispersion and strong Coulomb interactions. |37] 

We choose the microscopic parameters of the model Hamiltonian (|5|) such that a small 
number of experimentally measurable properties are fitted. These are the plasma frequencies 
and the Fermi surface shapes. We then proceed to calculate the spin fluctuation spectra and 
compare with experimental results. 
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B. Renormalized Band Structure 



The band structure of the Coulomb renormalized quasiparticles is described by the large 
N mean field theory. The mass enhancement is determined microscopically by the renor- 
malization of the hybridization matrix element, V P d — > V* d . The renormalization of the 
hybridization in mean field theory is mainly determined by the parameters associated with 
the copper states: the bare hybridization V p d and the bare level difference e p — e d (U is taken 
as infinite). We choose these bare parameters such that the renormalized plasma frequency, 
shown in Fig. |3](a), fit the experimental values given in Fig. ^(b) (derived from Drude fitting 
the optical conductivity spectra fl38|). 

The second important feature of the renormalized band structure is the shape of the Fermi 



surface. From the angular resolved photoemission spectroscopy (ARPES) measurement [39 



and the LDA band structure calculations j4(J, the Fermi surface in YBaCuO is found to 
be rotated by 45° relative to that of a diamond shape expected in a nearest neighbor tight 
binding form. The Fermi surface in LaSrCuO retains a diamond-like shape, though it 
is somewhat "twisted" so that the nearly flat regions of the Fermi surface are closer to 



the T point. |4(| The precise mechanism for the difference in the Fermi surfaces in these 
systems is at present not known. Since the spin degrees of freedom in the large U limit are 
mainly associated with the copper states, within our model study we choose appropriate 
oxygen dispersion to derive the corresponding shape of the Fermi surface for each family. 
Specifically, the Fermi surface "twisting" in LaSrCuO is derived by choosing a nonzero t pp , 
and the Fermi surface rotation in YBaCuO by choosing a non-zero t' pp in addition to t pp . 
These Fermi surfaces are shown in Fig. || 

The third feature associated with the renormalized band structure, specifically related to 
the two dimensionality, is the logarithmically divergent van Hove singularity in the density 
of states. The separation between the energy at which the van Hove singularity occurs and 
the Fermi level, u>vh = \Evh — Ep\, will be shown to be manifested in the spin fluctua- 
tion spectra. In this regard, we emphasize two important features. First, strong Coulomb 



10 



interactions pin the van Hove singularity near the Fermi level over a wide range of dop- 



ing concentrations. pT| , |42| Secondly, due to the different shapes of the Fermi surfaces in 
YBaCuO and LaSrCuO, the energy scale ujvh for YBaCuO is larger than that of LaSrCuO. 
Typically, we found ujvh of the order 25meV in YBaCuO, and less than 5meV in LaSrCuO. 



C. Magnetic Interactions and Dynamical Spin Susceptibility 

In the strong Coulomb correlation limit, the quasiparticles are dominantly associated 
with the copper states. Therefore, the exchange interaction between the quasiparticles can 
be thought of as primarily between the copper states, mediated by high energy states mainly 
of oxygen character. Such a picture represents a smooth evolution to doped case from the 
superexchange interaction in the half-filled limit, and can be formally established through 
analyzing the fluctuations beyond the mean field theory. ||TT| Our analysis shows that the 
overall amplitude of the magnetic interaction can be determined unambiguously only in the 
limit of large or vanishing t pp . Because t pp in the cuprates is in the intermediate region, we 
have appealed to phenomenology to pin down the amplitude of the magnetic interactions. 
Our general conclusion will be that, in the metallic region the magnetic interaction is of 
moderate strength. It is non-zero, but far from causing a magnetic instability. 

The q-dependence of the magnetic interaction can be specified on general grounds. We 
find that, for the antiferro magnetic interactions the q-dependence is to a good approximation 
given by Eq. Such a q-dependence reflects the smooth evolution from the superexchange 
interaction in the half-filled case. It is derived, as in the half-filled case, from the absence of 
the on-site term due to the exclusion of double occupancy as a result of the strong on-site 
Coulomb correlations. We emphasize that, these arguments refer to the character of high 
energy states, and do not depend on the details of the low energy quasiparticle dispersion. 
Indeed Eq. qualitatively describe the q-structure of the magnetic interaction in both 
YBaCuO and LaSrCuO systems, despite their very different Fermi surface shapes. 

The renormalized band structure determines the bare dynamical spin susceptibility as- 
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sociated with the quasiparticles— the Lindhard function. 

( \ ^ /(k)-/(k + q) 

where E(k) is the quasiparticle energy dispersion. In our approach, the Lindhard function 
has already incorporated substantial amount of interaction effects, through the Coulomb 
renormalization of the quasiparticle dispersion. It incorporates low energy scales induced by 
interaction effects, in addition to those associated with specific bandstructure such as van 
Hove singularity and nesting. The energy scale associated with the mass enhancement is 
called the coherence energy T co h. 

For a magnetic interaction which is not too close to induce a magnetic instability, the 



dynamical spin susceptibility has the generalized RPA form, [43 



x(q,<*>) ~ i — tt~\ — i — V 8 

i - ^(q)Xo(q,^) 

In this way, the magnetic interaction induces further softening of the antiferromagnetic 
fluctuations. When the magnetic interaction approaches to leading to a magnetic instability, 
processes beyond RPA will become important. This will be further discussed in Section |V|. 



IV. THEORY OF SPIN DYNAMICS IN THE METALLIC CUPRATES: 
NEUTRON SCATTERING AND NMR/NQR RELAXATION 

We now apply this formalism to the analysis of the spin fluctuation spectra. We aim at 
understanding the qualitative aspects of the experimental data. 

A. Wave Vector Dependence of the Dynamical Susceptibility in LaSrCuO 

In order to understand the neutron cross section, we start with a discussion of the 
Lindhard function in two dimensions. As is emphasized in Ref. 1T2], the (dynamical) Kohn 



anomaly in two dimensions leads to peaks in x"(q, uS) at the wave vectors q = 2kp, in 
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contrast to three dimensional case. This feature was emphasized independently in Ref. [15 . 



It has been further studied in Ref. [[0J within the one band Hubbard model for weak to 
intermediate Coulomb correlations. In an isotropic system, the Kohn anomaly gives rise to a 
ring of peaks. For a general Fermi surface, the shape of the peaked region in a; ) reflects 
the geometry of the Fermi surface. When applied to LaSrCuO with a Fermi surface given in 
Fig. H](a), nesting-like effects further enhance some portion of this peaked region. In fact, it 
induces four absolute maxima, at Q* = ^(1±<5, 1), -(1, 1±5). (|5| This is clearly seen in the 
three-dimensional q-structure shown in Fig. f|(a). The value of the incommensurability, 5, 
reflects the nesting wave vector which depends on the relative distances between the nearly 
flat regions of the Fermi surface: 5 increases when the nearly flat portions of the Fermi 
surface are closer to the T point. 

The magnetic interaction has a q-dependence of the form Eq. (^) and is peaked at 
Qaf = (- > -)• The effect of the magnetic interaction is to enhance 5(q, uS) around the whole 
region near Qaf- For moderate values of J a j J c , (where J c is the strength of the interaction 
which gives rise to a magnetic instability), the four peak structure is enhanced without 
changing shape. This is clearly seen in the three-dimensional q-structure at J Q /J C = 0.6 
shown in Fig. |^(b). When projected along two of the four peaks, as is shown in Fig. ^|(c), 
this can be directly compared to the experimental result shown in Fig. 0(b) (when resolution 
broadening is taken into consideration). 

We now turn to the implications of the experimentally found large and strongly dop- 
ing (x) dependent incommensurability (5) shown in Fig. 0(c). We plot in Fig. |7| 5 vs. 
x calculated within the present three-band large U model. For comparison, we also show 
the corresponding results for the small U Hubbard model with nearest neighbor hopping t 
only, and with next nearest neighbor hopping t' in addition to t. For the small U one band 
Hubbard model with only a nearest neighbor hopping, the incommensurability vanishes at 
half-filling and is strongly doping dependent: the nesting wave vector is equal to Qaf at half- 
filling, and deviates from Qaf as doping is increased. However, the incommensurability is 
smaller than the experimental value by more than a factor of two. When an additional next 
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nearest neighbor hopping is included, the Fermi surface is "twisted". With appropriately 
chosen t', the flat regions on the Fermi surface are further apart, and the incommensura- 
bility is enhanced. However, the next nearest neighbor hopping in this case destroys the 
special perfect nesting feature of the half-filled limit, and doping induces only very weak 
modifications to the incommensurability. This is again inconsistent with experimental 
results. 

In the three band large U case, the incommensurability is enhanced for the same reason 
that the flat regions of the Fermi surface are further apart compared to the one band nearest 
neighbor tight binding case. However, this change of the Fermi surface shape occurs due to 
the oxygen dispersion (the oxygen-oxygen hopping t pp ), which couples to the copper states 
through the renormalized hybridization V* d . In the strong U limit, the rigid band picture 
breaks down: the renormalized hybridization V* d decreases as the doping concentration 
is decreased. As a result, the amount of the Fermi surface "twisting" gets smaller, and 
the incommensurability is decreased. In the asymptotic half-filling limit, the renormalized 
hybridization V* d is reduced to zero, and the antiferromagnetic wave vector approaches the 
commensurate limit Qaf- Therefore, our results for the incommensurability is consistent 
with both features of the experimental data for incommensurability, i. e. the relative large 
values and the strong doping dependence. 

We emphasize that, in our calculation of 5*(q, uj) at different dopings, the bare param- 
eters of the Hamiltonian are fixed. The change in the renormalized band structure as the 
doping concentration changes is determined self-consistency from the large U constraint. 
Therefore, we argue that the experimentally observed large and strongly doping dependent 
incommensurate structure is a manifestation of the strong Coulomb correlation effect in the 
spin fluctuation spectra in LaSrCuO. 
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B. Wave Vector Dependence of the Dynamical Susceptibility in YBaCuO 



We now proceed to study the spin dynamics in the YBaCuO family. The q-structure of 
the Lindhard susceptibility corresponding to the appropriate band structure with a Fermi 
surface of Fig. |5](b) is shown in Fig. §(a). Because of the Fermi surface rotation, the Kohn 
anomaly induced— and nesting enhanced— peaks are far away from the antiferromagnetic 
wave vector Qaf = (-,-)• Therefore, the dynamical susceptibility is essentially featureless 
near Qaf- 

Since the magnetic interaction is peaked around Qaf, it enhances the amplitude of 
S^q, uj) or Jmx(q, uj) in the region surrounding Qaf- This results in an essentially commen- 
surate peak shown in Fig. ||(b). When projected along the diagonal direction, shown in Fig. 
H(c), our results can again be directly compared to experimental results such as shown in 
Fig. 0(a). Our assumption that J is not particularly close to leading to a magnetic insta- 
bility yields a peak relatively broad with a width essentially independent of temperature, as 
is illustrated in Fig. |(c). 

An exchange interaction peaked at antiferromagnetic wave vectors is therefore crucial 
to our analysis. Such a q-dependence for the magnetic interaction occurs naturally in the 
strong coupling limit, as emphasized in last section. This is very different from the weak 
coupling limit: the dominant exchange interaction is simply the on-site (q-independent) 
term U for small U. 

We conclude that, while the broad commensurate peaks of S(q,u) in YBaCuO are in 
striking contrast with the sharp incommensurate structure in LaSrCuO, they are in fact a 
consistency check within our theory: a) the absence of the incommensurate peaks in YBaCuO 
reflects the fermiology of this system, as does the presence of the incommensurate peaks in 
LaSrCuO; b) The presence of the commensurate peak in YBaCuO reflects an exchange 
interaction peaking at Qaf- This reflects the strong Coulomb correlations, as does the large 
and strongly doping dependent incommensurability in LaSrCuO. 

A number of other cuprate families have rotated Fermi surfaces. Our analysis leads to 
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the prediction that spin fluctuations in these systems will also be commensurate. 



C. Temperature and Frequency Dependence of the Dynamical Spin Susceptibility 

The role of low energy scales has been the theme of many studies in the metallic cuprates. 
|T0| Low energy scales can arise from features in the bandstructure such as the van Hove 
singularity {ojvh defined above) and/or the proximity to nesting. They can also occur due 
to the proximity to various instabilities; we call the soft energy scale due to the proximity 
to the Mott localization the coherent energy (T co h) and that due to the proximity to the 
magnetic instability the spin fluctuation energy (w s f). 

The frequency and temperature dependence of the dynamical susceptibility in the present 
theory for both YBaCuO and LaSrCuO are discussed extensively in Ref. [14]]. The lowest 



energy scale is manifested as a weak peak in the frequency dependence, and is found to 
correspond to the van Hove singularity energy ujvh defined above. This energy scale is larger 
for YBaCuO than that for LaSrCuO. At a moderate strength, the magnetic interactions 
slightly lower the peak energy, and considerably enhance the low energy spin fluctuation 
spectra. 

Alternatively, the low energy scales can be illustrated in a scaling plot of the susceptibility 
in terms of uj/T. When a low energy scale, say E*, is present, the frequency dependence of 
the spin fluctuation spectra depend on uj in terms ofu/E* for frequencies and temperatures 
smaller than E*. On the other hand, at higher frequencies and temperatures, the spin 
fluctuation spectra can scale with uj/T. A signature for a low energy scale is a deviation 
from Lo/T scaling at low temperature and frequencies. In Fig. || such scaling plots are 
shown for susceptibilities calculated for both LaSrCuO and YBaCuO systems. A deviation 
from scaling is clearly seen in YBaCuO, when the frequency is smaller than ujvh and the 
temperature smaller than Tvh ~ \^vh, where ojvh ~ 25meV. This prediction is consistent 
with subsequent experimental results of Ref. |]25|](a) for YBaCuO shown in Fig. |3|. The 



scaling behavior persists to lower energies in LaSrCuO, since ujvh is smaller. [I7[] Systematic 
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experimental studies in LaSrCuO is at this stage under way. 
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D. NMR/NQR Relaxation Rate 

Thus far, we have given a systematic presentation of the wave vector, frequency, and 
temperature dependence of x"(q, cu), as is reflected in neutron scattering data for both 
YBaCuO and LaSrCuO. We now turn to the analysis of the NMR/NQR data. 

Shown in Fig. |lO|(a) is the calculated temperature dependence of the NMR relaxation rate 
at the copper site in YBaCuO case, which can be directly compared with the experimental 
data above T c given in Fig. |I[ The results for LaSrCuO are similar. The high temperature 
saturation in our calculation occurs due to a proximity to Mott localization and the moderate 
magnetic interaction effects. The proximity to Mott localization is reflected in the moderate 
mass enhancement, leading to a coherent energy scale T co h above which the spin fluctuation 
spectral weight starts to saturate. The antiferromagnetic interaction leads to a further 
softening of the spin fluctuation spectra. The van Hove energy scale discussed in the context 
of neutron data plays a relatively minor role in 1/Ti due to the integrability of the logarithmic 
divergence. 

Since considerable interaction effects are already incorporated in Xoi °ur assumption of a 
moderate strength of the magnetic interaction J is consistent with both the experimentally 
observed enhancement and the deviation from linear in temperature dependence in 
This is in contrast with weak coupling approaches in which the interaction has to be fine- 
tuned to be extremely close to leading to a magnetic instability in order to yield soft-enough 
temperature scale in the temperature dependence of 57- and strong enough enhancement in 
its magnitude. 

The situation for {jr)o is more complex. Because peaks in the the dynamical spin sus- 
ceptibility is strongly incommensurate in LaSrCuO, and nearly-commensurate but broad in 
YBaCuO, a perfect form-factor cancellation is not expected in either system. Indeed, as can 
be seen in Fig. [IT](b) for the YBaCuO case, (^r)o shows considerable deviation from the 
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linear in temperature dependence. We note that, experimental study of (^-)o is at present 
less thorough than that for (tjt-)cu- Recent data indicate that, over certain tempera- 
ture range, (ijr-)o i n nearly optimal doped LaSrCuO may not be inconsistent with what is 
expected from an incommensurate spin fluctuation spectra as revealed in the neutron scat- 
tering. The reconciliation of neutron and (^-)o deduced spin fluctuation spectra is among 
the most important issues to be further studied, both experimentally and theoretically. 



V. FURTHER DISCUSSION 

The spin dynamics in the cuprates is an extensive subject. We have given a detailed 
description of certain aspects. In this section, we sketch on some related issues. 

A. Additional Aspects of the Spin Dynamics in the Normal State 

In this subsection, we briefly comment on the double layer magnetic coupling in YBaCuO, 
the temperature dependence of the uniform spin susceptibility, and the related "pseudo-gap" 
effects in the dynamical spin spectra. 

One important question in the physics of the high T c copper oxides is the role of interlayer 
couplings in the low energy electronic dynamics. In this regard, one feature which has 



been established experimentally is the double-layer magnetic coupling in YBaCuO. [49| It 
is possible to address this problem microscopically in terms of Coulomb renormalizations. 
As a first step towards understanding the phenomenological implications, we have carried 
out an RPA analysis of this double-layer magnetic correlation in Ref. ]n|(b). Incorporating 
an interp layer electronic coupling (t±) extracted from the splitting in the ARPES-induced 
Fermi surfaces |HJ in YBa2Cu30y_5, we find that both the frequency dependence and the 
c-direction wave vector q z dependence of S(q,u) can be understood provided an interlayer 
magnetic coupling (J±) is present. The extracted J± and t± are comparable in magnitude, 
suggesting that non-perturbative effects can be important. 
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Some classes of cuprates exhibit a gap-like features in the magnetic dynamics above T c . 
This is most extensively studied in the deoxygenated YBaCuO in which the spin susceptibil- 
ity x as well as (5^)0 an d (^t)y (and over a much narrower temperature range (^r^)c u ) all 
decrease substantially as temperature is decreased, 0] and a gap-like feature occurs in the 
frequency dependence of the dynamical spin susceptibility [pq] . Along the line of approach 
taken here, higher order processes beyond RPA become important when the magnetic in- 
teraction drives the system close to an antiferromagnetic instability. When these additional 
processes are taken into account, a new energy scale can be generated below which x de- 



creases as temperature is lowered. |5T] In addition, specific quasiparticle dispersion can lead 
to a "pseudo-gap" in S(q,u) for more extended q values. |52[ Detailed analysis is needed 
to see whether the "pseudo-gap" features indicate additional non-perturbative effects. A 
discussion of these "pseudo-gap" features in terms of RVB singlet formation can be found in 
Ref. p3[. Recently, it has been proposed that these "pseudo-gap" features may be associated 



with charge degrees of freedom which occur for all q, and are related to the short coherence 
length in the superconducting states. []54| Alternative arguments have been made that, they 
are associated with the interlayer singlet formation between adjacent CuC>2 layers within a 
unit cell. []5S| 



B. Spin versus Charge Dynamics in the Normal State 

One central issue in the metallic cuprates is the relation between spin and charge dynam- 
ics. This is highlighted by the proposal of spin charge separation. |56| Experimentally, the 
low energy spin fluctuation spectra exhibit various low energy scales. As we have discussed in 
detail, they are well described by the coherent spectra associated with the correlated quasi- 
particles. In addition, the spin fluctuation spectra exhibit crossover behavior at intermediate 
energy range (which we loosely define as, say, a decade of scale around lOOmeV). These 
features are strikingly similar with spin fluctuation spectra in the heavy fermion systems. 
[ TD| , |5T]] In charge responses, however, low energy scales appear only in non-optimally-doped 
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cuprates while intermediate energy crossover does not seem to occur, as demonstrated in 
the ubiquitous linear temperature dependence in resistivity and, perhaps more strikingly, 
the quadratic temperature dependence in the Hall angle over a wide range of temperatures 
|59| . At higher energies, not much is known about spin dynamics in the metallic cuprates, 
while various anomalous features in the charge dynamics have been found through optical, 
Raman, and photoemission spectroscopy studies. A systematic study comparing spin and 
charge dynamics at various energy scales should reveal much on the physics of the copper 
oxides. 

In the half-filled Mott insulator, strong Coulomb correlations push the spectral weight for 
charge excitations to higher energies, leaving spin excitations to play a dominant role at low 
energies. It is therefore natural to address how the low energy spin fluctuations affect the low 
energy charge dynamics in the doped case. Several groups have addressed the contribution 



of the spin fluctuations to resistivity. p7| , f41| In the context of neutron data, we mentioned 
that the lowest energy scale in the spin fluctuation spectra within our model came from the 
two-dimensional van Hove singularity, wyn- in Ref. [M, we found that ujvh is small over 



a wide range of doping concentration due to the strong Coulomb correlations. This leads 
to a linear in T dependence in the spin fluctuation induced scattering rate, over a wide 
temperature range above Tyu- Furthermore, since ujvh plays only a minor role in the NMR 
relaxation rate due to the integrability of the logarithmic divergence, this is consistent with 
differences in the low energy scales observed in different experimental probes. Further work 
is needed to address whether the spin fluctuations alone can account for all the anomalous 



low energy charge responses, or additional anomalous charge excitations [§^,^3] have to be 
invoked. 

C. Spin Dynamics in the Superconducting State 

Considerable attention has focussed recently on the spin dynamics below T c . The central 
issue has been the extent to which one can infer, from the spin fluctuation spectra, the nature 
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of the superconducting order parameter in the system: anisotropic versus isotropic pairing 
and the existence or not of nodes at the Fermi surface. Theoretical interpretations of the 
Cu and O Knight shifts and NMR relaxation rates in YE^CusOy when combined with 
anisotropy studies |Kj have been argued to give strong support for a d x 2_ y 2 pairing state. 
pT| , |52] ] In addition, recent neutron data on both Lai.8 5 Sr 15 Cu0 4 and YBa2Cu 3 66 reveal 
temperature T and frequency u dependent features which show considerable spin fluctuation 
spectral weight within the superconducting gap. p3 j^4|,p5 | 

Taking into consideration the interplay between fermiology and the Coulomb correlation 



effects, we have found that |63| the anomalous temperature dependences at low frequencies 
observed in neutron measurements p3,24|,25| of S(q, u) are compatible with a d x 2_ y 2 pairing 



state in both LaSrCuO and YBaCuO. The calculated change in the q-structure of S(q, uj) as 



one goes from above to below T c , due to a d x 2_ y 2 order parameter p2| , is expected to occur 
most strikingly in YBaCuO. The lack of such a change in q-structure in recent experiments 



in both YBaCuO |25| and LaSrCuO []64| poses a serious question to the anisotropic pairing 



scenario. 



VI. CONCLUSIONS 

The spin fluctuation spectra in the normal state of the metallic cuprates exhibit various 
anomalies. The contrasting wave vector dependence of the dynamical spin structure factor 
for YBaCuO and LaSrCuO, as well as their anomalous temperature and frequency depen- 
dence, provide the basis for our argument that, both the fermiology and strong Coulomb 
correlations play important roles in the low energy spin dynamics. 

We have shown that, the existence in LaSrCuO and the absence in YBaCuO of incom- 
mensurate antiferromagnetic peaks are natural consequences of the different shapes of the 
Fermi surfaces. At the same time, both the large value and the strong doping dependence 
of the incommensurability in LaSrCuO, and the presence of the commensurate peaks in 
YBaCuO, are the manifestation of the strong Coulomb correlation effects (along with fermi- 
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ology effects). Within our analysis, the strong Coulomb correlations lead to two effects. It 
renormalizes the quasiparticles, and induces residual exchange interactions. Because of the 
quasiparticle renormalization, we find that there can be enough spin fluctuation spectral 
weight to explain the temperature and frequency dependence in the neutron and NMR re- 
laxation rate, with an antiferromagnetic interaction which is far from being on the verge of 
leading to a magnetic instability. Meanwhile, the various energy scales found in the resulting 
spin fluctuation spectral function are consistent with energy scales within the renormalized 
coherent quasiparticle description of the spin dynamics. 

We are thus led to conclude that, the low energy spin dynamics of the metallic cuprates 
are appropriately described in terms of correlated quasiparticles with a Luttinger Fermi 
surface, in contrast to excitations in a frustrated spin system. These quasiparticles interact 
with a moderate value of residual antiferromagnetic interaction. These conclusions place 
strong constraints on the issue of the relevance of the antiferromagnetic spin fluctuations to 
superconductivity. They also provide a prototype for the low energy spin dynamics in 
the doped Mott insulators. 
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FIGURES 

FIG. 1. Spin lattice relaxation rates a) at the planar copper and oxygen sites in YBa2Cu307_,5 



(with T c = 93K) [18] and b) at the planar copper site in Lai.85Sro.isCu04 [17|. 



FIG. 2. Neutron scattering cross section as a function of momentum q a) in YBa2Cu306.6 at 
T=10K for a diagonal scan q = h^-(l, 1) (The q dependence for T > T c is essentially the same. 
From Ref. |25|](a)); and b) in Lai.86Sro.i4Cu04 at T=35K for a scan q = \{Q X + ^,Qx — f), i-e. 
along two of the four maxima at Q* = f (1 ± 5, 1), f (1, 1 ± 5) where 5 = 0.245 (from Ref. @(a)). 
Plotted in c) is the incommensurability versus the doping concentration (from Ref. [p3|(b)). 

FIG. 3. The q-integrated dynamical spin susceptibility deduced from neutron scattering cross 
section in YBa2Cu30g.6 as a function of T/u. The arrows show the temperature below which 
x"^) deviates from T/lo scaling. Data of Ref. |25|](a). 

FIG. 4. Plasma frequency cv p as a function of hole concentration x a) calculated from the 



renormalized quasiparticle band [14] and b) derived from Drude-fitting the optical spectra [38]. 



The "fitted" curve in b) represents u> p ~ ^fx. 

FIG. 5. The Fermi surface of the renormalized quasiparticle band for a) Lai g5Sro.i5Cu04 and 
b) YBa2Cu307 (from Ref. [0]). The LDA-calculated Fermi surface is shown for c) Lai.sSro.2Cu04 
(from Ref. |^](a)) and d) YBa2Cu307 (shaded areas, Ref. |I(|(b)). Open circles in d) correspond 
to the ARPES-measured Fermi surface for YBa2Cu306.g (from Ref. |39|](b)). 



FIG. 6. Calculated S^q, uj) versus (q x ,Qy) for Lai.g2Sro.i8Cu04 a) from the Lindhard function 
(J = 0) contribution; and b) for J / J c = 0.6; Also shown is c) the projection of b) along two 
maxima, with q = ^ (k + | , k — | ) and 5 = 0.34. Here the temperature and frequency are 1 and 
10 meV respectively. From Ref. ]l2|(b). 



FIG. 7. Comparison of the incommensurability 5 as a function of doping concentration x in 
one band models (dashed line, ti = and dot-dashed line ^Ai = 0.25) with the present theory 
(solid circles). From Ref. p4|]. 
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FIG. 8. Calculated 5(q, u) versus (q x ,qy) for YBa2Cu3C"6.7 a) from the Lindhard function 
(J = 0) contribution; b) for J /J c = 0.7. Here the temperature and frequency are 1 and 10 meV 
respectively. Also shown is c) the projection of b) along the diagonal direction. From Ref. |l3|| (a). 

FIG. 9. Normalized x"(q, u) as a function of ui/T for a) Lai^Sro.isCuCM and b) YBa2Cu30y. 
The solid lines represent the scaling curves. From Ref. ]l3[(b). 



FIG. 10. Temperature dependence of NMR relaxation rate in YBa2Cu307_<5 (a) at the copper 



site and b) at the oxygen site. From Ref. [14 
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